Cosmo logical Dynamics from Modified f(R) 
Gravity in Einstein Frame 

Daris Samart 1 
Fundamental Physics & Cosmology Research Unit 
The Tah Poe Academia Institute (TPTP), Department of Physics 
Naresuan University, Phitsanulok 65000, Thailand 

Abstract 

In this paper, we investigate and analyze the cosmological dynamics of the 
universe, with an effect of modified f(R) gravity emerging at cosmological scale. 
We choose the Einstein frame as a physical frame. We consider phase portraits 
of the universe at the late time from modified f(R) gravity model. This re- 
sult gives our universe an acceleration phase expansion without introducing 
existence of dark energy dominating our universe. 
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1 Introduction 

One of the greatest problems in cosmology today is the problem of accelerated ex- 
pansion of the universe. Fortunately, we can directly obtain various data parameters 
of the universe from observations. The highly precise observational data from both 
high presifted supernovae type la [1] and the anisotropy power spectrum of CMB [2] 
exactly tell us that our universe is increasing its rate of expansion to more than that in 
the past. Almost every cosmologist believes that the universe must be dominated by 
some mysterious field called dark energy (for interesting recent review see [3]), which 
drives the accelerated expansion of the universe. An alternative explanation that 
solves this problem modifies Einstein's gravity theory (i.e. general relativity, or GR 
for short) to account for the accelerated expansion (for general review see [4]). Carroll 
et al. [5] proposed a modified model of gravity by adding an inverse power of Ricci 
scalar into the Einstein-Hilbert action. It resulted in a number of deviations from GR 
at small curvature and cosmological scale, in which the universe is observed to be ex- 
panding with acceleration. The model gives us a de Sitter- and an anti-de Sitter-space 
solutions in vacuum case, which provide us with a purely gravitation-driven universe, 
whose expansion is accelerated. It is therefore in good agreement with current ob- 
servational data. This model is approached by standard variational principle [5], i.e. 
the action is varied with respect to the metric only. The result of which is the fourth 
order differential equation. Unfortunately, the model suffers instabilities problem [6] , 
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and has been ruled out, despite its equivalence to scalar-tensor theory [7], by a solar 
system experiment considering a gravitational bound state. The problem seemed to 
be resolved after Vollick [8] used Palatini variational approach that treated the metic 
and affine connection as independent variables in order to derive the field equation, 
which is a second order differential equation [11] ; and later, Nojiri and Odintsov [29] 
solved the same problem by using metric formalisms in which R 2 is added into the 
action of the 1/R model. However, Ref. [12] commented that this model should have 
a conflict in itself, because it was also suitable for a particle experiment at small 
energy scale. But Ref. [13] pointed out that the method which was discussed by Ref. 
[12] was just a mathematical equivalence, not a physical one (for extended discussion 
see [10, 11]). After all, the interesting of this model is strengthened by the fact that 
it can be derived from string/M-theory, as shown by Ref. [14]. Several authors have 
worked on various aspects of the model and extended it [16, 21, 22, 29]. However, 
almost all of these works used the model to find modified versions of Friedmann 
equation [8, 21, 24, 29] to fit with or be constrained by such observational data as 
SNa Type la, CMB anisotropy. Only a handful of works [18, 19] did investigate or 
analyze the dynamical system of this model. Ref. [18] considered the phase portraits 
of H vs H (time derivative of Hubble parameter and itself) in Jordan frame, whereas 
[19] considered between R vs i?(time derivative of Ricci scalar and itself). We will 
consider cosmological dynamics of modified gravity with variables identical to that of 
Ref. [19], but the metric variation formalism in phase portrait of Einstein frame is 
used instead. In our scenario, the Einstein frame is the physical frame, which gives 
us self-gravity of the scalar field's effective potential V(<p) (i.e. <fi in Einstein frame 
representative form of Ricci curvature in Jordan frame) [17, 23, 26]. This paper is or- 
ganized as follows: the conformal transformation between Jordan frame and Einstein 
frame is briefly reviewed in section 2, plus the use of Einstein frame to consider f(R) 
model of gravity, in which the modified Friedmann equation is found by following the 
process of [30]; next the f(R) gravity model is investigated and analyzed in section 
3 by using cosmological variables obtained from the Friedmann equation. Finally we 
conclude with problems and the future direction for this paper. 

2 Conformal Transformation Approach to the Ein- 
stein Frame and Friedman Equation 

In this section, we review the conformal transformation method following [24, 25, 26] ; 
see [20] for the standard derivation. The action of f(R) gravity in Jordan frame with 
matter field can be written in the form 




(1) 
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where k 2 = 8irG, and Sm is the action of matter field. Ricci scalar is defined by 
R = g ab R a b , Rab = R c a cbi an d the curvature tensor is defined by 

rid a p f)~rd J_ F e F ' F e F ^ /o\ 

11 abc — u c L ab u b L ac ~t~ 1 at 1 ce 1 ac 1 6e ■ l z J 

Varying this action with respect to the metric tensor g ab , we get 
5S, 



Sg 



ab 



f(R)'R ab - -gabf(R) - V„V 6 /'(i2) + g ab V c V c f{R) = K 2 T ah , (3) 



where = df(R)/dR, and T ab = 5SM/8g ab is the energy-momentum tensor. 

Note that the field equation derived from (3) in Jordan frame gives us a fourth 
order differential equation, which is too complicated to solve [5, 18]. Magnano and 
Sokolowski [26] was argued the f(R) theories in vacuum case. There are is the 
existence of the Einstein frame that (we will define it below) is minimally coupling 
to the metric of the physical frame. We will analyze the late time universe where we 
assume universe is vacuum. So we choose the Einstein frame as a physical frame. The 
metric formalisms seem to be equivalent to Palatini formalisms, i.e. g a b is compatible 
with V a , therefore an afline connection is necessarily equal to the Christoffel symbol 
in general relativity (f(R) = R) [20]. But both variation methods are obviously 
different in the nonlinear function of the Ricci scalar (f(R) is a higher term which 
metric formalisms give rise to a fourth order equation but Palatini formalisms give 
rise to a second order equation [9, 11] ), and we obtain the conformal transformation 
between two frames as 2 

E 9 b = e+g ab , (4) 

where we introduce as a new scalar field defined by = In f'(R) [15, 23]. We can 
rewrite the action in (1), which is dynamically equivalent to the Jordan-Helmholtz 
action without matter field [14, 26, 27, 28], as 



^ J {f(v)+x{R-*)) V^d'x, (5) 

where x an d cr are auxiliary fields. After varying (5) with respect to x, we get 

a = R, (6) 
and after varying (5) with respect to a again, 

X = /V) • (7) 

X is eliminated by using (7), hence 



SjH '- 2^ 



l ~2 J {f(a) + f'(a)(R-a)) ^gd'x. (8) 



2 It's convenient to obtain the conformal factor by using Palatini formalisms where the action 
is varied with respect to T c ab . 
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Note that (8) is in the form of the scalar-tensor gravitational action for f(R) gravity 
(1) in Jordan frame [28] . We can rewrite the action (1) in Einstein frame by applying 
conformal transformation to (4), (i.e. g ab — > e^g ab ) [14, 26], 



S E [ E g ab ] = ^- 2 J ( E R-l V b V a 0V 6 0-y(0)) V^9d 4 x, 



(9) 



where V a is the covariant derivative compatible with E g a b, and V(0) is the effective 
potential which depends on a choice of f(R) models. The effective potential is defined 

= Rf>(R) - f(R) 

mr ( ' 

Varying the action (9) with respect to E g ab , we get 

E R ab - \ E g ab E R = 3V a V 5 - \ E g ab Q E g ab V a <t> V fe + , (11) 

whilst varying it with respect to , we obtian 

V^0+y = O, (12) 

where = dV/dcfi. Now the present universe is considered based on FRW metric, 
which, in spatially flat space, has the form 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) , (13) 

where a(t) is a scale factor of the universe. 
Now we will calculate the modified Friedmann equation in Einstein frame from the 
line element (13) and the (0,0) component of the field equation (11). The Friedmann 
equation then becomes 

where H = a/a is the Hubble parameter, a = da/dt, and (15) is used to obtain (14). 
We also directly obtain the equation of motion for the scalar field from (12) in 
Einstein frame metric, i.e. 

+ 3#0+^ = 0. (15) 

In the next section, we will consider cosmological dynamics in phase portraits of some 
particular types of f(R) gravity model. 



3 Our effective potential form is equivalence to the Palatini formulation [17] . 
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3 Cosmological Dynamics of f(R) Gravity 



In this section, we will consider cosmological evolution in the phase portrait of Einstein 
frame, with self gravitating scalar field, to investigate and analyze the acceleration 
phase for some types of f(R) gravity models. We follow the analysis of [30]. 
Ref. [29, 31] was demonstrated mechanism of inflation start/end and acceleration 
start via the potential of curvature evolution. Therefore, we can analogue this mech- 
anism by manifest the scalar field dominates the late time universe in the Einstein 
frame (the curvature is very small). The acceleration can happen if the field at late 
time is moving very slowly (i.e. the curvature decrease very slow). Therefore, we can 
use the slow-roll approximation [33] in this phase [30], i.e. 0^0 and <j) 2 <C V . We 
use the slow-roll approximation only to find the late time trajectory in the phase 
space, both the Friedman equation and and the equation of motion of the scalar field 
are 

H 2 * ^, (16) 
3H<P ~ (17) 

By applying the slow-roll approximation condition to both the Friedmann equation 
in (16), and the equation of motion for the scalar field in (17), then we directly get 
the late time trajectory in the phase space portrait as 

0~-^L. (18) 
3V6V 

Now we derive the accelerating condition of the universe in f(R) gravity model. The 
accelerating condition of the expansion is 

- = H + H 2 >0, (19) 
a 

which implies that 

-^<!" (20) 

We obtained H directly by differentiating H 2 in (16), and together with using (17), 
we then obtain 

1 v 2 

if = - — (21) 
18 V K J 

By using (16) and (21) substituting in (20) , it's becomes 



V 2 
v 



< 3 . (22) 
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These equations are the accelerating conditions of expansion. In the previous section, 
we have derived the Friedmann equation in Einstein frame metric. We will derive 
the dynamical equation of <fi and . For obtain new value of H with out slow-roll 
approximation by using (14) and (15) 

H = ~l<P 2 - (23) 

We define a new variables as 



X = <j>, Y = <j>. (24) 

Using new variables in (24) substitute into (23) and (15) , we obtain autonomous 
system as 

H(X,Y) = - 3 -Y* 
X = Y 

Y = - m x,Y)Y-l d -^f>. (25) 
The function H(X, Y) is given by 



H(X,Y) = J\y* + ±V(X). (26) 

These equations will be used for the analysis in the phase portrait of the f(R) modified 
gravity. Next, we will investigate and analyze the phase space portraits for some types 
of the effective potential V(<f>), which depends on the choice of f(R) models of gravity. 
By using the autonomous system of equations in (25) . We also show the acceleration 
condition curve which divides acceleration and non-acceleration regions in the phase 
space portraits. 

3.1 f(R) = R- n 2 ( n+ V /R n model 

This model was proposed by [5] (for n — 1) at the cosmological scale when curvature 
is very small at the late-time universe. Our gravitational Lagrangian density (n — 1) 
in Jordan frame has 

f(R) = R-^/R, (27) 

where fi 4 is a positive-valued parameter. The effective potential V(<p) in Einstein 
frame depends on the choice of f(R) ■ At the late time, R is very small, we get 

e* = m) = l + £*£. (28) 




Figure 1: The phase portrait of <j) vs (f) in the modified 1/R gravity model when /j, = 1 . 
We note that everywhere in phase space portrait is an acceleration region following the 
condition in (32) . 

Then, the effective potential in the Einstein frame of the model f(R) = R — p r/R 
gravity at the late time by using in (28) substituted into (10) in [29] 

K($~A e -§*. (29) 
By plugging the effective potential (29) into (18) , then the late time trajectory is 

' e-5*. (30) 



2v"3p 

The effective potential of our late time approximation looks like the power law infla- 
tion, which gives scale factor evolution proportional to the power of time, a(t) ~ t p . 

This proportionality comes from the potential 4 V ~ exp ^— | j with p = | in our 

case [33, 34] , such that 

a ~ t 4/3 . (31) 

We note that potential increases when curvature is minimum and therefore the uni- 
verse begins to have the power law expansion and substitute (29) into (22), then 
accelerating condition is 



d I 2 ' ' 




< 3 . (32) 



4 Our convention compared with standard power law formalisms is tfi = \pn t f> with potential taking 
the form V ~ exp(— w |yj) . 
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Figure 2: The phase portrait of <fi vs <f) in the modified lni? gravity model, when a = 1 , the 
acceleration region divided by acceleration condition (blue line) where acceleration region 
is bright region and shaded region is no acceleration region. 

The above condition obviously shows that everywhere in the phase space portrait of 
vs in f(R) = R — fi 4 /R gravity model can be an acceleration region, i.e. wherever 
any attractor occurs in phase space, it can always generate the accelerated expansion 
in both present and late time phase and the phase plot will occur attractor along the 
line of the late time trajectory curve (30) ;See Figure 1. 

3.2 f(R) = R + a In (R/'y 2 ) model 

This model was proposed by [31] to solve the problem of instabilities in the fourth 
order differential equation from CDTT's model in metric formalisms [5] . It also 
provides accelerating phase at late time. (For more detail of Palatini formalisms of 
this model, see [32]). The gravitational Lagrangian density of this model in Jordan 
frame is 

f(R) = R + alniR/^ 2 ), (33) 

where a and 7 are positive-valued parameters. As in the previous subsection, at the 
late time R is very small, we have 

= f(R) = 1 + I ~ I (34) 

then we obtain the effective potential in the Einstein frame of the model f(R) = 
R + aln(i?/7 2 ) gravity at the late time by plugging (34) into (10), as in [31] 

V(0)~ae- 2 ^0. (35) 



S 



The late time trajectory is 



■ ge-'(l - 20) 

0= 3^ • (36) 



This effective potential of our late time approximations that looks like generalized 
exponential potential was originally examined by J. D. Barrow [35], and extended 
investigation was done by [36] , i.e. the potential takes the form V ~ m exp(— A0 n ) . 
For our case, m — n — 1 which gives scale factor evolution that is proportional to the 
power of time as [31] 

a(t) ~ ti . (37) 

We performed calculations following the previous subsection by plugging (35) into 
(22), the accelerating condition is then 

2 - 20 - 1 < . (38) 

The above condition manifests itself as a curve (blue line) dividing the acceleration 
region (bright region) from the no acceleration region (shaded region) . We note that 
from Figure 2, an attractor occurs in the phase plot in acceleration region, therefore, 
the universe's expansion will be accelerated in f(R) = R + a In (-R/7 2 ) gravity model, 
and the attractor occurs in the phase plot along the line of the late time trajectory 
curve in (36) ; See Figure 2. 



4 Conclusions 

We have investigated the modified f(R) gravity models, namely 1/R and lni? gravity 
models. Both models give the late time attractors. These results mean that our uni- 
verse has its expansion accelerated. The evolution when scalar field is at some given 
point on the potential has to be independent of any initial conditions. The attractor 
of 1/R gravity can drive the universe accelerating its expansion, because the condition 
of acceleration allows everywhere in phase space portrait to be an acceleration region. 
In the lni? gravity, an attractor occurs in the acceleration region according to the 
acceleration condition. The universe can then accelerate its expansion in this model. 
It is surprising that our analysis in Einstein frame gives similar result as in [37] which 
chose Jordan frame as a physical frame. They showed that f(R) model gravity has a 
late time attractor on a phase space portrait of H vs H. Both Jordan and Einstein 
frames are analyzed similarly for f(R) gravity model. Our investigation and analysis 
give well behavior in cosmological dynamics at late time. The f(R) gravity theories 
also provide inflation at early time by adding higher curvature in the gravitational 
action [29, 38] . These results are confirmed by observational data in several works 
that we have discussed above. It also passes the solar system experimental tests, 
and gives no stabilities problem [29, 31, 39] . Hopefully, the modified f(R) gravity 
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is a viable candidate to dark energy. Although our analysis in Einstein frame gives 
late time acceleration phase, Carroll et al.[5] originally analyzed universe in Jordan 
frame. They considered evolution in Einstein frame and transformed it back to Jor- 
dan frame as the physical frame. Because, in Einstein frame, the energy density p 
evolution with scale factor p ~ a~ 3 does not hold in the ordinary matter dominated 
case 5 , which is the problem that will be investigated in the future work. However, 
we cannot tell whether the Jordan frame or Einstein frame is a physical frame. For 
general consideration and mathematical equivalence (via conformal transformations) 
between Eintein, Jordan and even ideal fluid frames of f(R) gravity was given in 
detail in [10] . 

Note Added During the preparation of this work. Modified gravity which seems 
to be quite successful as dark energy was criticized by recent papers from L. Amen- 
dola, D. Polarski and S. Tsujikawa appeared on arXiv [40] because some simplest 
version of it like 1/R n gravity are not able to describe at once and matter dominated 
phase, then transition from matter dominated phase to acceleration, then acceler- 
ation. However, some counterexamples where such sequence of universe epochs is 
possible in more complicated class of f(R) gravity ware suggested in S. Capozziello, 
S. Nojiri, S. D. Odintsov and A. Troisi [41] . 
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